Lower bounds on the complexity of simulating quantum gates 
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We give a simple proof of a formula for the minimal time required to simulate a two-qubit unitary 
operation using a fixed two-qubit Hamiltonian together with fast local unitaries. We also note that 
a related lower bound holds for arbitrary n-qubit gates. 



I. INTRODUCTION 

Understanding quantum dynamics is at the heart of 
quantum physics. Recent ideas from quantum compu- 
tation have stimulated interest in studying the physical 
resources needed to implement quantum operations. In 
addition to a qualitative understanding of what resources 
are necessary, we would like to quantify the resource re- 
quirements for universal quantum computation and other 
information processing tasks. Ultimately, we would like 
to understand the minimal resources that are necessary 
and sufficient to implement particular quantum dynam- 
ics. 

As a first step towards answering these questions, it 
has been shown that there is a sense in which all entan- 
gling dynamics are qualitatively equivalent. In particu- 
lar, it has been shown that any n-qudit two-body Hamil- 
tonian capable of creating entanglement between any pair 
of qudits is, in principle, universal for quantum compu- 
tation, when assisted by arbitrary single-qudit unitaries 
[HElESEIiSli- Thus > any particular entangling 
two-qudit Hamiltonian can be used to simulate any other, 
provided local unitaries are available. This suggests that 
such dynamics are a fungible physical resource. 

Having established the qualitative equivalence of all 
entangling dynamics, we would like to quantify their in- 
formation processing power. In particular, it is interest- 
ing to consider the minimal time required to implement 
a unitary operation, U, on a two-qubit system, using a 
fixed Hamiltonian, H, and the ability to intersperse fast 
local unitary operations on the two qubits. This problem 
was studied by Khaneja, Brockett and Glaser |£g, who 
found a solution using the theory of Lie groups. Their 
results, although giving a solution in principle, are nei- 
ther explicit about the form of the minimal time, nor do 
they explain how to construct all elements of the time- 
optimal simulation. Further work by Vidal, Hammerer, 
and Cirac 10], from a different point of view, resulted 
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in an explicit formula for the minimal time, and gave 
a constructive procedure for minimizing that time (see 
also where an alternate proof is given by the same 
authors). 

The purpose of the present paper is to give a simplified 
proof that the formula of Vidal, Hammerer and Cirac 
is, in fact, a lower bound on the simulation time. Note 
that the difficult part of [Tol Hl| was proving the lower 
bound; finding a protocol to meet the lower bound was 
comparatively easy. 

The main advantages of our proof are its simplicity 
and conceptual clarity, as compared to the ingenious, but 
rather complex, arguments in [iJElO]. This simplicity 
is achieved by making use of a powerful result from linear 
algebra, Thompson's theorem. We expect that Thomp- 
son's theorem might be useful for many other problems 
in quantum information theory. A second advantage of 
using Thompson's theorem is that it does not rely on 
special properties of two-qubit unitary operators. There- 
fore, essentially the same arguments give a lower bound 
on the time required to implement an n-qubit unitary 
operation using a fixed n-qubit interaction Hamiltonian, 
and fast local unitary operations. 

Our approach to the proof of the lower bound has its 
roots in the framework of dynamic strength measures for 
quantum operations |12(. The dynamic strength frame- 
work is an attempt to develop a quantitative theory of the 
power of dynamical operations for information process- 
ing. The idea is to associate with a quantum dynamical 
operation, such as a unitary operation U, a quantita- 
tive measure of its "strength." In it was shown that 
such strength measures can be used to analyze the min- 
imal time required for the implementation of a quantum 
operation. The present paper takes a similar approach, 
but instead of using a single real number to quantify dy- 
namic strength, we use a vector-valued measure. This 
can also be compared to the analysis of optimal sim- 
ulation of Hamiltonian dynamics using a set of several 
strength measures [13|. 

Our paper is structured as follows. Section ITT1 reviews 
some background material on majorization, Thompson's 
theorem, and the structure of the two-qubit unitary ma- 
trices. The main result, the lower bound on optimal 
simulation, is proved in Section II I II We conclude in 
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Section IIVI by presenting our generalization of the lower 
bound to n qubits and suggesting some directions for fu- 
ture work. In addition, an appendix gives a procedure 
for calculating a canonical decomposition of two-qubit 
unitary gates. 

II. BACKGROUND 

This section reviews the relevant background needed 
for our proof. Section III Al reviews the basic notions of 
majorization, introduces Thompson's theorem, and ex- 
plains how to use Thompson's theorem and majorization 
to relate properties of a product of unitary operators to 
properties of the individual unitaries. Section III Bl in- 
troduces the canonical decomposition, a useful represen- 
tation theorem for two-qubit unitary operators, and Sec- 
tion lll Ol presents an analogous decomposition for Hamil- 
tonians. 



A. Majorization and Thompson's theorem 

Our analysis uses the theory of majorization together 
with Thompson's theorem. More detailed introductions 
to majorization may be found in |14|. Chapter 2 and 3 
of O, and in [H[l3. 

Suppose x = [x\, . . . , xd) and y = (j/i, . . . , Vd) are two 
.D-dimensional real vectors. The relation x is majorized 
by y, written x -< y, is intended to capture the intuitive 
notion that x is less ordered (i.e., more disordered) than 
y. To make the formal definition we introduce the nota- 
tion { to denote the components of a vector rearranged 
into non- increasing order, so x^- = (x{, . . . ,x\j), where 
x'^. Then x is majorized by y, that is, 

x -< y, if 

i=i 3=1 

for k = 1, . . . , D — 1 , and the inequality holds with equal- 
ity when k = D. 

To connect majorization to Hamiltonian simulation, 
we use a result of Thompson relating a product of two 
unitary operators to the individual unitary operators. 
Recall that an arbitrary pair of unitary operators can be 
written in the form e lH and e lK , for some Hermitian H 
and K. Thompson's theorem provides a representation 
for the product e lH e lK in terms of H and K: 

Theorem 1 (Thompson [18]). Let H,K be Hermitian 
matrices. Then there exist unitary matrices U, V such 
that 

e iH e iK = e i(UHU'*-+VKV'<)_ 

The proof of Thompson's theorem in [131 depends on 
a result conjectured earlier by Horn A proof of 



this conjecture had been announced and outlined by Lid- 
skii [20| at the time of Thompson's paper. However, re- 
marks in 0| suggest that [2(j did not contain enough 
detail to be considered a fully rigorous proof. Fortu- 
nately, a proof of Horn's conjecture has recently been 
fully completed and published. See, for example, |2lLl22| 
for reviews and references. 

Thompson's theorem may be related to majorization 
using the following theorem of Ky Fan: 

Theorem 2 (Ky Fan OHl). Let H, K be Hermitian 
matrices. Then \{H + K) -< X(H) + X(K), where X(A) 
denotes the vector whose entries are the eigenvalues of the 
Hermitian matrix A, arranged into non-increasing order. 

Combining the results of Ky Fan and Thompson, we 
have the following: 

Corollary 3. Let H, K be Hermitian matrices. Then 
there exists a Hermitian matrix L such that 

e lH e lK = e lL - \{L) -< X(H) + X(K). (3) 

We will not apply this corollary directly, but we have in- 
cluded it here because it captures the spirit of our later 
argument, combining the Thompson and Ky Fan theo- 
rems to relate the properties of a product of unitaries 
to the individual unitaries themselves. Corollary [3] can 
be regarded as a vector-valued analogue of the chaining 
property for dynamic strength measures used in [l2T | to 
establish lower bounds on computational complexity. 



B. The canonical decomposition of a two-qubit gate 

The canonical decomposition is a useful representa- 
tion theorem characterizing the non-local properties of 
a two-qubit unitary operator. It was proved by Khaneja, 
Brockett, and Glaser using ideas from Lie theory. 
Kraus and Cirac |24| have given a constructive proof 
using elementary notions, while Zhang et al. |25j have 
discussed the decomposition in detail from the point of 
view of Lie theory. The decomposition states that any 
two-qubit unitary U may be written in the form 

U={A X ® B 1 )e* e " x * x+e v Y9Y+9 ' z ® z )(A 2 ® B 2 ), (4) 

where A\, A2, B±, B2 are single-qubit unitaries, and the 
three parameters, 9 X , 9 y , and 8 Z characterize the non- 
local properties of U. 1 Without loss of generality, we 
may choose the local unitaries to ensure that 

~>e x >e y >\e z \, (5) 



Prior to 9], Makhlin 26 1 gave a proof that the non-local proper- 
ties of U arc completely characterized by 8 X , 8 y and 8 Z , but did 
not write down the canonical decomposition explicitly. 
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and we refer to the set of parameters chosen in this way as 
the canonical parameters for U . We will see below that 
these parameters are unique. We define the canonical 
form of U to be 



(A\®Bl)U(At®Bl) 



(6) 



up to local unitaries, U c is equivalent to U. It will be 
convenient to assume through the remainder of this sec- 
tion that U has unit determinant. This is equivalent to 
requiring that A\, Ai,B\,Bi can all be chosen to have 
unit determinant. 

The canonical parameters turn out to be crucial to 
results about simulation of two-qubit gates. If 



U _ e i(8 x X®X+8 v Y(g>Y+9 z Z®Z) 



(7) 



is the canonical form of U, then we define the non-local 
content, <j>(U), of U by <f>(U) := X(Hjj), where 

H v := 9 X X ® X + 6 y Y ® Y + 6 Z Z ® Z. (8) 

Explicitly, the components of 4>(U) are 



(9) 
(10) 

(11) 
(12) 



We now outline a simple procedure to determine the 
canonical parameters of a two- qubi t unitary operator. 
Our explanation initially follows [Tlj and [23. However, 
as explained below, there is an ambiguity in the proce- 
dure described in those papers, related to the fact that 
the logarithm function has many branches. Our proce- 
dure resolves this ambiguity. 

To explain the procedure, we need to introduce a 
piece of notation, and explain a simple observation about 
single-qubit unitary matrices. The spin flip operation on 
an arbitrary two-qubit operator is defined as 



M 



(Y ® Y)M T (Y ® Y), 



(13) 



where Y is the Pauli sigma y matrix, and the transpose 
operation is taken with respect to the computational ba- 
sis. Note that the spin flip operation may also be written 
as M = M T , where the transpose is taken with respect 
to a different basis, the magic basis |28j . 



|QQ) + |ii) 

■ |01) + |10) 



• 100) 



111) 



V2 
101} - 110) 



V2 



(14) 



The observation about single-qubit unitary matrices that 
we need is the following. Let U be any single-qubit uni- 
tary matrix with unit determinant. Then 



uyu t = y, 



(15) 



where the transpose is taken in the computational basis. 
This simple identity is easily verified. 

Now suppose U is an arbitrary two-qubit unitary with 
unit determinant. By definition of the spin flip, and sub- 
stituting the canonical decomposition, we have 

UU = (A 1 <S>B 1 )U C (A 2 ®B 2 )(Y<S>Y) 

•x(Al®B 2 T )U c {AT®Bl)(Y®Y). (16) 

By the identity Eq. I|15|) we see that 

UU = (Ax <» Bx)U c (Y <g> Y)U c {Al ® Bf)(Y <g> Y). (17) 

Using the fact that Y ®Y commutes with X<8)X, Y®Y, 
and Z ® Z, we see that Y (8 Y commutes with U c , and 
thus 

UU = [At ® B x )Ul{Y ® Y)(AJ ® Bf )(Y ® Y). (18) 

Finally, applying Eq. (fT^) again gives 

UU = {A 1 ®B 1 )Ul{A\®B\) (19) 

Eq. I|19|l suggests a procedure to determine the canonical 
parameters for U, based on the observation that 



X(UU) = A(C/ C 2 ) = (e 2 *,e 2, *,e 2 '> 3 ,e 2#4 ), 



(20) 



where the <pj are related to the canonical parameters 
6 X , y and 8 Z by Eqs. lj§)l~ H12[) . It is tempting to conclude 
that one can determine 9 x ,9 y , 9 Z from the eigenvalues of 
UU, simply by taking logarithms and inverting the re- 
sulting linear equations. Indeed, such a conclusion is 
reached in and [53 , using arguments similar to those 
just described. Unfortunately, determining the canon- 
ical parameters is not quite as simple as this, because 
z — > e lz is not a uniquely invertible function. In partic- 
ular, e lz = e l ( z + 2 ^ m ) ; where m is any integer, so there 
is some ambiguity about which branch of the logarithm 
function to use in calculating the canonical parameters. 
In fact, we prove later that no one branch of the logarithm 
function can be used. However, these considerations do 
allow us to reach the following conclusion: 



Lemma 4. Let U be a two-qubit unitary, 
exists a Hermitian H such that 



Then the 



UU 



X(H) = 0(C/). 



Moreover, if H is any Hermitian matrix such 
X(UU) = \{e 2iH ) then it follows that X(H) = <j>{U)-{ 
where rh is some vector of integers. 



(21) 

that 
■Km, 



Although this lemma is sufficient to prove our later re- 
sults, there is in fact a simple method for exactly calcu- 
lating the canonical parameters. Because there are many 
applications of the canonical decomposition, we describe 
this method in the appendix. The method will not be 
needed elsewhere in the paper. 
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C. The canonical form of a two-qubit Hamiltonian 

Finally, we introduce one additional concept, the 
canonical form of a two-qubit Hamiltonian, H [3j . Any 
two-qubit Hamiltonian H can be expanded as 

3 

H = h jk aj®a k . (22) 
j,k=o 

Then let 

H -4- // 

H' := = J2 h i k a o ® ° k - ( 23 ) 

That is, H' is just the Hamiltonian that results when 
the local terms in H are removed. It is not difficult to 
show that H and H' are interchangeable resources for 
simulation in the sense that, given fast local unitaries, 
evolution according to H for a time t can be simulated 
by evolution according to H c for a time t, and vice versa. 
Furthermore, by doing appropriate local unitaries, it can 
be shown [3| that simulating H' (and thus H) is equiva- 
lent to simulating the canonical form of H, 

H c = h x X®X + h y Y®Y + h z Z®Z, (24) 

where h x > h y > \h z \. Once again, H and H c are inter- 
changeable resources for simulation. 

Note that the three parameters h x ,h y ,h z are com- 
pletely characterized by the three degrees of freedom in 
A(if c ) = X(H + H)/2, just as the three canonical param- 
eters x ,8 y , 9 Z are completely characterized by the three 
degrees of freedom in X(U 2 ) = X(UU). 

III. SIMULATION OF TWO-QUBIT GATES 

We now return to the main purpose of the paper, prov- 
ing results about the time to simulate a unitary gate us- 
ing entangling Hamiltonians and fast local gates. We aim 
to prove the following result: 

Theorem 5 (Vidal, Hammerer, Cirac PI EH, cf. 
Khaneja, Brockett and Glaser [9]). Let U be a two- 
qubit unitary operator, and let H be a two-qubit entan- 
gling Hamiltonian. Then the minimal time required to 
simulate U using H and fast local unitaries is the mini- 
mal value of t such that there exists a vector of integers 
rh satisfying 

<j>(U) + Trm -4 A(g 2 + g) t. (25) 

Note further that only two vectors of integers need to be 
checked, rh = (0, 0, 0, 0) and rh — (1, 1, —1, —1), since all 
the other possibilities give rise to weaker constraints on 
the minimal time, t \wi Hi) . The difficult part of the 
proof of Theorem [S] is the proof that Eq. (|25|l is a lower 
bound on the simulation time, t, and it is this part of the 



proof that we focus on simplifying. The proof that this 
lower bound may be achieved follows from standard re- 
sults on majorization, and we refer the interested reader 
to [H 113 for details. 

To prove that Eq. I|25|) constrains the minimal time 
for simulation, we begin by characterizing the canoni- 
cal decomposition of a product of unitary matrices. Let 
A(U) := X(UU), and define the equivalence relation 
A ~ B for Hermitian matrices A and B iff X(A) = X(B). 
Then we have: 

Lemma 6. Let Uj be unitary matrices, and let Hj be 
Hermitian matrices such that Uj Uj — e 2 lHj . Then there 
exist Hermitian matrices Kj such that Hj ~ Kj, and 

A{U N . . . f/i) = x(e 2l(Kl+ - +KN) ). (26) 

Proof. We induct on N. The result is trivial for N = 1, 
so we need only consider the inductive step. Using the 
fact X(AB) = X(BA), we have 

A(U N+1 ...U 1 ) = X{U N+1 U N+1 U N ... UxUx . . . U N ). 

(27) 

By the inductive hypothesis there exist Hermitian Kj 
such that Hj ~ Kj and 

X(U N . . . Ux&i ...U N ) = X(e M{K 'i + - +K '^). (28) 

Therefore, U N . . . U x Ui . . . U N = e 2i ( K " + - +K '^ , for 
some Kj ~ Hj. Observe also that 

U N+1 U N+1 ~ U N+1 U N+1 - e 2iH »+\ (29) 

and thus Un+iUn+i = e 2tK '"+ 1 for some K'^ +l ~ H N+1 . 
It follows by substitution that 

A(U N+1 . . . Ux) = A(e 2lK «+ie 2l(K "+- +A « ) ). (30) 

Applying Thompson's theorem gives 

A{U N+1 ...U x ) = x(e 2 ^ Kl+ - +KN ^) (31) 

for some Kj ~ Kj ~ Hj, which completes the inductive 
step of the proof. □ 

Given this result, it is straightforward to complete the 
proof of Eq. 

Proof. Write U in the form 

U = e~ lHtl V ie ~ lHt2 V2 . . . V k - 1 e- lHtk , (32) 

where t±, . . . , t k are times of evolution, t = t± + . . . + i/. 
is the total time for simulation, and Vj are local uni- 
taries. Without loss of generality, we may assume H is 
in canonical form. Applying Lemma El we obtain 

A(U) = x{e 2l{H ^+- +Hktk) ) (33) 

where Hj ~ H for each j. Here we have used the ob- 
servation VjVj = 1, so all the contributions from local 
unitaries vanish. It follows from Lemma 0] that 

(t>(U)+nm = X(H 1 t 1 + ... + H m t m ), (34) 
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and using Ky Fan's theorem gives 

4>(U)+7rm^X(H)(ti + ... + t m ), (35) 
which is Eq. 1)25(1 . as desired. □ 



IV. DISCUSSION 

In this paper, we have provided a simplified proof of 
a lower bound on the time required to simulate a two- 
qubit unitary gate using a given two-qubit interaction 
Hamiltonian and local unitaries. The bound follows eas- 
ily from standard results on majorization together with 
Thompson's theorem on products of unitary operators. 

Although we have described canonical decompositions 
of two-qubit gates in some detail, we note that our proof 
does not actually require properties of the decomposition 
unique to two qubits. In fact, it is straightforward to 
prove an analogue of Eq. (12511 for an n-qubit system. For 
an n-qubit operator M, suppose we define a generalized 
spin flip M — > M, where M is the transpose operation in 
a basis such that, whenever M is local, M is orthogonal, 
i.e., MM — I. It is not difficult to construct examples of 
such bases, at least when n is even. An example is the 
basis obtained by rotating the computational basis us- 
ing the transformation (/ — iY® n )/y/2, for n even. This 
basis change gives M = Y® n M T Y® n , where the trans- 
pose is taken in the computational basis, and thus this 
operation generalizes the transpose in the magic basis. 
In this general setting the following lower bound on the 
time required to implement an n-qubit gate holds: 

Corollary 7. Let U be an n-qubit unitary operator, and 
let H be an n-qubit Hamiltonian. Then the time required 
to simulate U using H and fast local unitaries satisfies 

X - arg X(UU) + Trm -< X( - H + H ^ t . (36) 
for some vector of integers m. 

The proof follows simply by taking the arguments of both 
sides of Eq. I|33l) and applying Ky Fan's theorem. All 
steps leading up to Eq. I|33|l remain valid for n-qubit sys- 
tems using the above definition of the generalized spin 
flip. 

Unfortunately, we have not found any interesting ex- 
amples with n > 2 for which Eq. I|36(l provides a non- 
trivial lower bound on the time required to implement 
some quantum gate. It would be interesting to construct 
cases where Eq. (|36|l (or some similar condition) does 
give a nontrivial constraint on multipartite gate simula- 
tion. One might imagine that such techniques could be 
used to prove circuit lower bounds on certain quantum 
computations, although it does not seem likely that such 
bounds would be especially strong, given the well-known 
difficulty of this problem. 
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APPENDIX: A METHOD FOR COMPUTING 
THE CANONICAL PARAMETERS OF A 
TWO-QUBIT UNITARY GATE 

In this appendix, we describe a method for computing 
the canonical parameters of a two-qubit unitary, based 
on the discussion in Section III Bl The key is to take 
logarithms in just the right way. From Eqs. © and Q- 
(|T2"|l . we see that 

y > 20i > 20 2 > 203 > 204 > -y. (37) 

It is not difficult to find examples where the first or last 
inequality is saturated, so no single fixed branch of the 
logarithm function can be used to determine the (f>j . One 
might hope instead that there exists a method for choos- 
ing a different branch for each particular U, so that the 
corresponding 20^ lie within that branch. However, even 
this is not possible in general. To understand this, note 
that 

20! - 204 = 4(6 X + 6 y ). (38) 

In cases where 6 X — y = tt/4, we have 20i — 204 = 2tt, 
in which case the values 20-,- do not lie in any one branch. 

We now show how to compute the 0j. The idea is 
that we can first take the argument of the eigenvalues in 
Eq. H2U|) over some fixed branch. Then we can systemat- 
ically determine which of the resulting values have been 
shifted by 2tt from the value 20j (due to an incorrect 
branch) and correct these values accordingly. 

Let Sj, j = 1, . . . , 4 be defined as follows: 

2Sj = arg(e 2 ^)- (39) 

That is, 2Sj are the arguments of the eigenvalues of UU, 
where we take the argument over the branch {—■§, 
so that the Sj are contained in the interval (— f,4f]- 
Considering the range of values that <f>j may take, from 
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Eq. (|37|l . and the particular branch we are using, it is 
clear that: 

S- = i^ +n if ^<~f (-40) 

3 1 4>j otherwise. v ' 

From Eqs. l|5|l- lJT2j l we have 

01 + 02 + 03 + 04 = 0. (41) 

Combining Eqs. (|40|l and l|41(l . we see that 

S 1 + S 2 + S 3 + S 4 = irn, (42) 

where n is the number of 4>j that are less than or equal 
to —j. Possible values for n are 0,1,2 and 3 (all four 
(f>j cannot simultaneously be < — since that would 
contradict Eq. (|41fl h Since the 4>j obey the ordering in 
Eq. (|37() . then the n values of 0j that are less than or 
equal to — j are 04, . . . , 04_ n +i, and the remaining 4 — 
n values greater than j are 0i, . . . , 04- n - Thus, using 
Eq. (|40|l . we see that the set of values Sj consist of n 
"shifted" 4>j values 

04 + 7T, . . . , 04-n+l + 7T, (43) 



and 4 — n "non-shifted" values of 4>j 



0i,...,04-n- (44) 

Furthermore, all of the shifted values in l|43|) are no less 
than any of the non-shifted values in (|44|l . This is shown 
by combining Eq. JSJ with Eq. Ij38(l . giving 0i — 04 < 7r, 
which when combined with Eq. (|37|l implies that <j)j < 
0fc + 7r for all j, fc, as required. Therefore, the largest 
n values of Sj are guaranteed to be the values in l|43Jl . 
Thus subtracting ir from the largest n values of , gives 
us 04, ... , 04_„+i, and the the remaining 4 — n values of 
Sj give us 0i, . . . ,0 4 -„. 

In summary, the nonlocal parameters 9 X , 9 y and 9 Z may 
be computed as follows. Find the arguments of the eigen- 
values of UU over the branch (— f f Call these values 
2Sj. Calculate n = (Si + S2 + S3 + S^/n. Replace the n 
largest values of Sj by those values minus ir. The result- 
ing values, when placed in nonincreasing order, are equal 
to (0i, 02, 03, 04). The parameters 9 Xl 9 y and 9 Z are then 
found by inverting Eqs. (|5|)- (fT2^l . 
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